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We prove a theorem that shows the degeneracy of many-body states for particles in a periodic 
lattice and under a uniform magnetic field depends on the total particle number and the flux filling 
ratio. Non-interacting fermions and weakly interacting bosons are given as two examples. For the 
latter case, the phenomenon can also be physically understood in terms of destructive quantum 
interference of multiple symmetry-related tunneling paths between classical energy minima, which 
is reminiscent of the spin-parity effect discovered in magnetic molecular clusters. We also show 
that the quantum ground state of a mesoscopic number of bosons in this system is not a simple 
mean- field state but a fragmented state even for very weak interactions. 



Recently, cold atoms in optical lattices subject to a 
large effective magnetic field has become a new direction 
in cold atom research. Interacting bosons in rotating op- 
tical lattices have been investigated experimentally pQ, 
where rotation plays a similar role to magnetic field for 
charged particles. A synthetic magnetic field for neutral 
atoms has also been successfully realized through engi- 
neering atom-light interaction, and this technique can be 
applied to optical lattices straightforwardly [2]. These 
experimental progresses and the rich physics in such a 
system have generated lots of theoretical interest [3HS]. 
Most of the work so far are mean-field studies or classical 
Monte Carlo simulations of mean-field states [3-5 , and 
some focus on strongly correlated quantum-Hall states 
in the strong interaction limit [6]. In this letter, (i) we 
shall first prove a general theorem on the degeneracy of 
many-body states in this system, and (ii) we shall present 
non-interacting fermions and weakly interacting bosons 
as two concrete examples of this theorem. In addition, 
(iii) we will show that in the case of finite number of 
bosons, the system exhibits an intriguing non-mean- field 
ground state, even in the regime of very weak interac- 
tions. 

We consider a two-dimensional optical lattices and a 
uniform effective magnetic field B along z direction (cre- 
ated by rotation or atom-light interaction). In this work, 
we focus on the properties of a uniform system, and do 
not consider the harmonic trapping potential. Under the 
Landau gauge, the single particle Hamiltonian is given 

by 

Ho = ~dl + ^~ {-ihd y - eBxf + V olj {x,y). (1) 
2m 2m 

Ri and R2 are the basis vectors of the optical lattice 
potential Vol(x,v). The Hamiltonian for a many-body 
system reads 

N 

H = J2Ho(n) + J2 v (n-ri)- (2) 

i=l i<j 



Let N denote the total number of particles, and v = (f)/<fio 
denote the flux filling ratio, where <fi = £?(Ri x R2) • z is 
the effective magnetic flux per plaquette and </>o = h/e is 
the flux quantum. 

Theorem on Degeneracy. For v = p/q where p and 
q are coprime numbers, q is a prime number, and N/q 
is not an integer, all many-body eigenstates are at least 
q-fold degenerate. 

To prove this theorem, we first introduce the mag- 
netic translation operator for particle j as Tj(r) = 
exp{zrllj/^}, where Uj X = —ihd Xj — eByj and Uj y = 
—ihd Vj . One can show that both Tj(Ri) and Tj(R 2 ) 
commute with #o( r j)j but 

f i (R 1 )T i (R 2 ) = e^fj^fjCRi). (3) 

Furthermore we define the magnetic translation operator 
for all N particles together of one basis vector as 

Both Tx and Ty commute with J^i<j ^( r i — r j) smce a U 
particles are translated together. Therefore, we have 

[f x ,H} = [f Y ,H}=0. (5) 

However, 

f x f Y = e i27TpN/q f Y fx. (6) 

We can choose a common eigenstate of H and Ty as 
\&o with fy*o = 7?*o and H^ = e^ . Defining = 
(fx) l * , (Z = l,...,g-1), then 

f Y ^i = e- 2viplN/q r]^i, (7) 
H^i = e^t (8) 

If p and q are coprime numbers, and q is a prime number, 
when N/q is not an integer, pNl with I = 0, 1, . . . , q — 1 



are different integers modulo q, therefore all have dif- 
ferent eigenvalues of Ty, i.e., they are orthogonal states 
with the same energy. 

This theorem adds one more example to the rare cases 
that an exact theorem can be proved for a many-body 
system. It holds for any pair- wise interaction V(r), for 
both bosonic and fermionic systems, and for all lattice ge- 
ometry such as rectangular, triangular or hexagonal. It 
can also be generalized to a three-dimensional cubic lat- 
tice straightforwardly, since both Tx and Ty commute 
with translations along the z-direction. This theorem 
imposes a strong constraint on all the many-body the- 
ory of this system, namely, any approximation scheme 
applied to the system must respect the theorem. Here- 
after, we shall give two concrete examples using the two 
dimensional square lattice (Ri = ax and R 2 = ay) , from 
which we hope to provide a physical picture for this phe- 
nomenon. 

The first example is non-interacting fermions. For 
v =p/q, the eigenstates of Ho satisfy the magnetic Bloch 
theorem and are characterized by two good quantum 
numbers k x C (—7r/(aq),7r/(aq)] and k y C (— 7r/a,7r/a] 
[7] . Due to the magnetic translation symmetry, each sin- 
gle particle state is g-fold degenerate [5J |8]. Therefore, 
the ground state can not be unique unless TV is a multiple 
of q. While if N modulo q is I, the degeneracy is at least 
q\/{l\{q — 1)1) , which is a multiple of q. 

The second example is weakly interacting bosons. Tak- 
ing v — 1/3 as an example, with a simplified effective 
model, we find that the ground state is non-degenerate 
if N is a multiple of three, while three-fold degenerate 
otherwise. This is reminiscent of the spin-parity effect 
discovered in the studies of spin coherence of magnetic 
molecular clusters [9], where distinctive behaviors are 
found in spin integer and half-interger systems. There, 
the phenomenon is interpreted in terms of the quan- 
tum interference of multiple symmetry-related tunneling 
paths connecting the degenerate classical states. By sim- 
ilar quantum interference argument, we can understand 
the degeneracy for v = 1 /3 and more general cases where 
no analogy of spin-parity effect has been discussed before. 

Effective Model: Since each single particle state is q- 
fold degenerate, there must be g-fold degenerate single 
particle ground states. For v = 1/3, three degener- 
ate single particle ground states ^1,2,3 ( r ) are magnetic 
Bloch states with (k x ,k y ) = (0,0), (0,27r/(3a)) and 
(0, — 27r/(3a)). Let 61,2,3 be boson operators for these 
three modes, then all the states b\ n % n2 b\ {N ~ ni ~ n2) |0) 
(ni = 0, . . . , iV, n 2 = 0, . . . , N — n\) are degenerate with- 
out interactions. Interactions will result in two effects: it 
will mix these (N + 1)(N + 2)/2 states, and also intro- 
duce quantum depletion to other single-particle excited 
states. In the regime that the interaction strength is 
much smaller than the magnetic band width, the first 
effect is dominant due to the boson enhancement factor 
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FIG. 1: (Color online) This is a real space contour plot for the 
amplitude of the condensate wave function (<£i(r) + ^2(1*) + 
e l2ir/ V 3 (r))/V3 (a) and (^(r) + p 2 (r) + e -^ /3 ^ 3 (r))/V3 
(b). The intersections of the vertical and the horizontal 
lines indicate the lattice sites, i.e. the potential minima of 
Vol (x,y). The dark area indicates the locations of the vor- 
tices. 



and the absence of single-particle energy cost. In prac- 
tice, the band width is usually a small fraction (~ 0.01) 
of /z 2 /(2ma 2 ), while the interaction ~ h 2 a s /(ma 3 ), hence 
it requires that the scattering length a s should be at 
least three order of magnitude smaller than a ~ 0.5 jam. 
Atoms like 39 K, 7 Li and 133 Cs, which either have very 
small background scattering length or have the scattering 
length that can be tuned across zero by a Feshbach res- 
onance, are particularly suitable for reaching this limit. 

Hence, we take the approximation that only these three 
modes are kept, then the single particle term becomes a 
constant. Taking V(r) = U5(r), an effective Hamiltonian 
is purely given by the interaction as 

ffeff = a(h\ +h\+ n\) + 4/3(nin 2 + n 2 n 3 + hihs) 

+[27(S t 1 &I& 2 S 3 + 6+6+6x63 + 6^6162) + h.c], (9) 

where a = £/J>i| 4 d 2 r, (3 = U J |^i| 2 |(^ 2 | 2 ^ 2 r and 
7 = U j (p^^PscPr. A 90° rotational symmetry leads 
to a — 2(3 = 27. This model contains all the four-boson 
interaction terms allowed by momentum conservation. 

Mean-field Analysis. We first revisit the mean-field 
states (MFS) in the weakly interacting regime. Previ- 
ous mean-field studies of this system are limited to the 
Gutzwiller mean- field ansatz in tight-binding limit [3]. 
Here we implement another MFS written in the magnetic 
Bloch state bases, which are 



(U!b\ 



u 2 bl 



u 3 bl 



N 



|o>, 



(10) 



Three of them are 

(61 



By minimizing the energy we find six degenerate MFS. 

j = a] N \0)/y/N\ (j = 1,2,3), where 
«i = l&i + 6 2 + e l27r ^ 3 bs)/V3, d 2 and a 3 follow by cyclic 
permutation of the coefficients in a±. The other three 
dt N \0)/y/N\ (i = 1,2,3), where di = (61 +60 



are &j 



, 3 |0)/VAH (j = 1, 2, 3), where d 1 = (61 -r o 2 - 
e~* 27r / 3 6 3 )/<\/3, d 2 and ds follow by cyclic permutation of 
the coefficients in d\. 
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FIG. 2: (a) AE = (E 1 -E ) /N and (E 2 -E )/N as a function 
of N. E is in unit of U. (b): the eigenvalues of Xi/N of the 
density matrix p of the ground state. 



The vortex lattice configurations are displayed in Fig. 
[I] Here we obtain the magnetic Bloch function (pi(r) 
by diagonalizing the lattice potential Vol m the lowest 
Landau level subspace [7 j. These results agree with the 
Gutzwiller mean-field studies [3] and the classical Monte 
Carlo simulation in the tight-binding limit [4 . As one 
can see from Fig. [I] the unit cell is enlarged to q x q with 
vortex lattices, and therefore the symmetry of the vortex 
lattice state is lower than that of the original Hamilto- 
nian. Hence, the MFS are degenerate. The degeneracy 
of the MFS is also a manifestation of our theorem in the 
thermodynamic limit, and this degeneracy is crucial for 
the later discussion of the quantum ground state. 

Exact Diagonalization Studies. Given the total number 
of bosons N, we diagonalize the Hamiltonian Eq. [9] using 
the Fock state bases. We sort the eigen-energies from 
the lowest to the highest as E n . For the ground state, we 
compute the density matrix p = (b\bj) and its eigenvalues 
A^. The main results from this exact diagonalization are 
presented in Fig. [2]and[3j 

In Fig. |2|a) we plot AE = E 1 - E and E 2 - E for 
different N, which shows that the ground state is three- 
fold degenerate for N = 3n + 1 and 3n + 2, while it is not 
degenerate for N = 3n. Nevertheless, it also shows that 
AE decreases exponentially as TV increases. As we will 
show later, AE is in fact the tunneling splitting due to 
macroscopic quantum tunneling between different MFS 
with the same classical energy. In the limit N — >• +oo the 
difference vanishes and all states are practically at least 
g-fold degenerate. This numerical calculation clearly ver- 
ifies our theorem. 

In addition, in the case of a non-degenerate ground 
state, we show in Fig. [2jb) that all three eigenvalues A^ 
equal to N/3. Hence, the non-degenerate ground state is 
a fragmented state, instead of a single condensate. Tak- 
ing v — 1/3 and N = 24 as an example, the struc- 
ture of the quantum ground state is further illustrated 
in Fig. [3] We first rotate the single-particle bases to 
{aj(j = 1, 2, 3)}. The wave function is written as 



y 1 ^ni,n 2 

m ,ri2 



1 2 3 

i/ni!n 2 !(A r - ri\ - n 2 )l 



|0> (11) 




o o 







FIG. 3: (a) |*4ni,n 2 | as a function of (711,712) (b) \T>m 1 ,m 2 \ as 
a function of (7711,7712). See the definition in the text. Here 
N — 24. 



N 


D 


Wave functions 


3n 


1 


^1 + ^2 + ^3 + e i7riV/6 ($i + $ 2 + $3) 


3n+ 1 

or 
3n + 2 


3 


e -i2.N/3y i + ^ 2 + ^ 3 + V3e- l7riV/2 $i, 
1 + e - t2lTN/3 ^2 + ^3 + V3e- t7TN/2 <$>2, 
^1 + ^2 + e- l2lTN/3 ^3 + VSe- l7TN/2 <S> 3 



TABLE I: For v — 1/3, the relation between N, the degener- 
acy (D) and the ground state wave functions. 



|*4ni,n 2 | is plotted in Fig. [3^ a), where the most pro- 
nounced feature is the three peaks located at (711,712) = 
(iV, 0), (0, iV) and (0,0). It means that the ground state 
contains an equal weight superposition of three MFS ^ j 
(j — 1,2,3). Another feature in Fig. |3ja) is that there 
is also weight of A ni ,n 2 distributed broadly in the center 
region in the plot. This part of wave function can be 
recast as 



E 

mi ,rri2 



d[ ml d 



^/mi\rri2l(N — mi — 777-2)! 



|0>. (12) 



\^m u m 2 1 is plotted in Fig. [3^b), which again shows three 
sharp peaks located at (7771,7772) = (A/ - , 0), (0, A") and 
(0,0). Fig. [3] indicates the quantum ground state of 
this case can be well approximated (with overlap > 90%) 
as a coherent superposition of six degenerate MFS with 
appropriate coefficients. For a general situation of v = 
1/3, the detailed structures of the ground state are found 
from numerical calculations and are illustrated in Table 
I. 

Such a fragmented state results from quantum tunnel- 
ing between degenerate classical MFS. During the past 
years, quite a few of such examples have been studied in 
cold atoms systems. Fragmentation has been studied for 
bosons in double- well [10 , multi-well [11], and also in a 
single trap [12] or in a triple- well with dipolar interactions 
[T3] , The exact ground state of spin-1 bosons in the "po- 
lar" phase is found to be a singlet-pair condensate [T4] . 
Breakdown of mean-field approximation and fragmented 
ground states have also been found in a rotating BEC 
at critical frequency of vortex nucleation [I5j [16] and in 
the fast rotating regime of bosons [TJ. These states, in 
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principle, can be distinguished from MFS through mea- 
surement of noise statistics and interference pattern [10] . 
Nevertheless, they are usually very fragile and is difficult 
to prepare, in particular, for large N when the tunneling 
splitting between different MFS becomes exponentially 
small. However, the decay rate of this tunneling split- 
ting depends on the details of each model. It is still 
possible that there exists a regime that the decay rate 
is sufficiently low and a fragmented state with a good 
number of particles can be achieved experimentally. We 
leave this for future studies. 

Quantum Interference Viewpoint of Degeneracy. Pre- 
vious studies of spin coherence in magnetic molecular 
clusters have revealed an effect known as the spin-parity 
effect [9 . Consider a SU(2) spin model with time- 
reversal symmetry, and two classical energy minima lo- 
cated at the south and north poles of the Bloch sphere. 
One can show that for each tunneling path between these 
two classical minima, there is always another path re- 
lated by time-reversal symmetry, and therefore these two 
paths have equal tunneling amplitudes. One can also 
show that there is a relative phase exp{z27r5'} between 
the two paths. Therefore, the total tunneling ampli- 
tude is proportional to 1 + exp{z27r6'}, which vanishes 
for half-integer S and results in double degeneracy [9] . A 
two-mode boson model can be mapped to a SU(2) spin 
model with S = N/2. Recently, a similar effect has also 
been discussed for BEC in double-well [18] and should 
also exist in the model discussed in Ref. [12 . By the 
Schwinger boson scheme, H e ^ of Eq. [9] is mapped to a 
SU(3) spin model. However, to the best of our knowl- 
edge, no such effect has been reported in a SU(N) spin 
model for N ^ 3 before. The discussion below is equiv- 
alent to a generalized spin-parity effect for a SU(3) spin 
model, and similar analysis can even be possibly applied 
to a general SU (N) case. 

Using the coherent state representation, we write 



N 



\ Zj ). 



(13) 



where \zj) = e 



-\ Zj \ 2 /2 p z 3 a^ 



e Zja j\0). The tunneling amplitude 



between and ^ 2 can be formulated in terms of the 
imaginary time coherent state path integral as [9 



U 
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(*i|e- 



-Hr/h 



vn x 



( - \N 
7T Z N\ 



where S[Q] is the Euclidean action for path Q from the 

initial coherent state \z±) to the final one \z 2 ). We note 

that Hi is invariant under two different magnetic trans- 
. . 3 
lation operations: 



(i) b^bse-* 2 ^ 3 , h 2 ^b x e^'\ b 3 ^b 2 ; 

(ii) b x ^b 2 e-^'\ b 2 ^b 3 , h^b^'\ 



J27T/3 



under which di — >• e ±i27r/3 ai, d 2 — >• e ±i27r/3 a 2 , and a 3 — >• 
as (upper signs for (i) and lower signs for (ii)). Hence, 
for any tunneling path Q° (connecting initial state \z®) 
and final state l^)), there must be two other symmetry- 
related paths Q 1 and Q 2 (connecting \z\) and \z l 2 ) (i = 
1,2), respectively), which satisfy S[tt 2 } = S^ 1 } = SKP] 



and 



1,2 



e ±i2n/3 z and z h2 = eT i2*/3 z m By Eq 
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■a*Wu° 2 , where 



one has U{ 2 = e i27TN / 3 U? 2 and U? 2 = e 
U{ 2 (i = 0, 1, 2) denotes the tunneling amplitudes of path 
fi\ Thus, the total tunneling amplitude is always propor- 
tional to 1 + e i27rAr / 3 + e - i27rAr / 3 . For N = Sn + 1, 3n + 2, 
this tunneling vanishes due to destructive interference. 
This is a generalization of 1 + e l27TS (S = N/2) for the 
spin-parity effect. Similarly, one can show the tunneling 
between all and Sfrj (i 7^ j), and the tunneling between 
&i and <3>j (i ^ j) exhibits the same behavior. Hence, we 
have established an alternative viewpoint from quantum 
interference as to why the case TV being a multiply of 
three is different. 
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